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Nomenclature 

C  Coupling  function,  Eq.  (2) 

F  Gauss's  hypergeometric  function,  Eq.  (3) 

H  Dynamic  transfer  function 

H  Normalized  dynamic  transfer  function 

Nr  Number  of  blades  on  rotor 

V  Steady  axial  velocity 

a0 

Qq  Steady  volume  flow  rate  through  the  pump 


i  -  1 

f  Frequency  of  fluctuations,  Hz 

n  Rotational  speed  of  turbopump 

t  Time 

u  Velocity  of  the  rotor  cascade 


Ap  Unsteady  pressure  rise  through  the  rotor 

K 

Ap  Unsteady  pressure  rise  through  the  stator 

AV  Amplitude  of  axial  velocity  fluctuations 

a 

p  Density  of  the  liquid 

X  Stagger  angle  of  rotor 

R 

o  Solidity  of  rotor  =  chord  of  blade/spacing  of  blades 

R 

a)  Reduced  frequency  ■  f  cos  X  /N  $  n 

R  R  R 

v  Frequency  of  fluctuations,  radians/sec 

*  Flow  coefficient  *  V  /n 

a0 

r  Gamma  function 
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Introduction 

Nonsteady  cavity  flow  tests  conducted  in  a  closed  circuit  water  tunnel 
are  thought  to  induce  pressure  surges  in  the  tunnel  circuit  thus  invalidating 
the  test  results  [2].  There  has  been  sufficient  interest  generated  [2,4,5]  to 
warrant  understanding  and  analyzing  the  water  tunnel  circuit  as  a  dynamic 
system.  Indeed,  although  many  investigators  have  observed  and  mentioned  the 
existence  of  pressure  surges  in  closed  circuit  water  tunnels  [2,4,5],  there 
appears  to  be  no  systematic  study  that  has  been  performed  to  understand  the 
possible  existence  of  pressure  surges  and  to  quantify  the  surge  phenomenon  (if 
it  exists)  so  a  correction  procedure  can  be  adopted  for  nonsteady  cavity  flow 
test  data  obtained  from  closed  circuit  water  tunnel  tests. 

The  components  of  the  water  tunnel  circuit  that  affect  the  dynamic 
behavior  of  the  tunnel  circuit  have  been  identified  [3]  as: 

(i)  Unsteady  cavity  motion  in  the  test  section. 

(ii)  Phase  shift  in  the  unsteady  disturbance  due  to  the 
tunnel  circuit. 

(iii)  The  dynamic  characteristics  of  the  tunnel  drive  pump. 

(iv)  The  pressure  control  loop  involving  the  differential 
pressure  transducer  at  the  test  section,  the  pressure 
regulating  tank  and  the  controllers. 

In  this  memorandum,  the  dynamic  characteristics  of  the  48-inch  water 
tunnel  drive  pump  are  determined  using  the  theory  of  unsteady  flow  around  a 
two-dimensional,  linear  cascade  of  airfoils  developed  by  Ohashi  [1].  The 
results  obtained  show  some  interesting  features  that  could  be  important 
in  the  dynamic  analysis  of  the  water  tunnel  circuit. 
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Description  of  the  Tunnel  Drive  Pump 

The  axial  flow  pump  of  the  48-inch  water  tunnel  at  the  Garfield  Thomas 
Water  Tunnel  consists  of  a  four  bladed  impeller  95  inches  in  diameter  with 
flow  straightening  vanes  ahead  of  and  behind  the  impeller.  The  impeller  is 
driven  by  a  variable  speed  motor.  The  speed  can  be  varied  continuously  from  0 
to  180  rpm.  The  blades  of  the  impeller  may  be  adjusted  over  approximately  28 
degrees  [6].  The  solidity  of  the  impeller  is  about  1.0. 

Analysis 

The  relationship  between  the  pressure  rise  across  the  pump  and  the  flow 
rate  is  indicated  by  the  characteristic  curve  of  the  pump.  This 
characteristic  curve  is  valid  only  when  the  pump  operates  in  the  steady  state 
condition,  that  is,  with  a  constant  flow  rate.  When  the  flow  rate  fluctuates, 
the  pressure  rise  will  not  respond  fast  enough  to  follow  along  the  steady 
state  characteristic  curve. 

The  dynamic  characteristics  describe  a  relation  between  the  fluctuating 
volume  flow  rate  at  the  inlet  to  the  pump  and  the  fluctuating  pressure  rise 
across  the  pump,  as  a  function  of  the  frequency  of  fluctuations.  Both 
amplitude  and  phase  relations  have  to  be  obtained  for  a  meaningful  use  of 
the  dynamic  characteristics. 

It  is  assumed  that  (1] 

(i)  The  flow  is  a  two-dimensional  cascade  flow  for  each  blade 
element  and  it  is  inviscid,  incompressible  and  non- 
cavitating. 

(ii)  The  angular  momentum  of  the  inlet  flow  into  the  rotor 
is  steady. 


(iii)  The  unsteadiness  is  periodic. 
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(iv)  The  increase  of  velocity  head  at  the  outlet  edge  of  the 
rotor  is  converted  into  pressure  rise  in  the  following 
stator  or  diffuser  without  any  time  delay. 

The  unsteady  pressure  rise  through  the  rotor  and  stator  is  given  by 

fl], 


%  +  Aps 

P 


u  AV 


tan  AR[l  -  exp(-  ttctr) ] 
1  +  iuiR  C(mR,oR) 


ivt 


where 

Ap  ■  Unsteady  pressure  rise  through  the  rotor 
K 

Ap  =  Unsteady  pressure  rise  through  the  stator 
S 

p  =  Density  of  the  liquid 

u  ■  Velocity  of  the  rotor  cascade 

AV  =  Amplitude  of  axial  velocity  fluctuations 
a 

A  -  Stagger  angle  of  rotor 
K 

o  =  Solidity  of  rotor  =  chord  of  blade/spacing  of  blades 
R 

ui  =  Reduced  frequency  ■  f  cos  A  /N  $  n 
R  R  R 

f  *  Frequency  of  fluctuations,  Hz 

v  »  2itf  *  Frequency  of  fluctuations,  radians/sec 

t  *  time 

i  -  -  1 

n  *  Rotational  speed  of  turbopump 

N  *  Number  of  blades  on  rotor 

R 

$  m  Flow  coefficient  -  V  /n 
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V  =  Steady  axial  velocity 
a0 

And  C(wR,aR)  the  coupling  function  given  by  [1]  is 


rl1  +  v  r(ia)R) 


c(v°r)  ’  -  k  t  *  .,r",  t  '  ‘1  +  :,i  .  -r  ,  W- 1  •  t  +  1ur  ; 


“r  r(i  +  in.)  r(i  +  iul 


exp(-  2ttor))  -  l]  ,  (2) 

where  T  =  Gamma  function  and  F(a,b,c,z)  is  the  Gauss's  hypergeometric  function 
given  by 

00  (a)  (b) 

F(a,b,c,z)  =  l  *—±zn  ,  c  *  0,-1, -2 .  (3) 

n=0  n 


where 


,  s  r (a  +  n) 
(a)n  =  “(a) 


n  »  T(b  +  n) 
(b  n  -  ~<b)  “ 


(C)  =  Hi 
^  ;n  rCc) 


The  dimensionless  form  of  unsteady  pressure  rise  is  given  by 


tpR  H,s  -  2  tan  XR[l  -  «p[-  «»R)]  i0  ivt 
l/2ou2  1  +  i«R  Ct»R,oRj  Qq 


-7- 


26  August  1985 
KR:lhz 


where 

AQq  =  Amplitude  of  unsteady  volume  flow  rate  through  the  pump 
Qq  =  Steady  volume  flow  rate  through  the  pump 

The  dynamic  transfer  function  H(ui)  is  therefore  given  by 

h„  .  „(t )  -  2  tan  XR[  1  -  exp(-  iroR)] 


H(u)  = 


R+S 


(AQ(t)/Q0j  1  +  itoR  C[a)R,oRJ 


(8) 


ivt 


where  AQ(t)  =  AQ^e  *  ^en  there  are  no  fluctuations  in  volume  flow  rate, 


the  quasi  steady  limit  is  given  by 


AQ, 


0 

(hR+s)  =  -  2  tan  Ap(l  -  exp(-  ircrp))^ —  $  exp(ivt)  . 


qS 


R  ;Qf 


(9) 


We  may  now  form  the  ratio, 


H(u>) 


R+S 


1 


^hR+S^  1  +  1U)R  C*“r’°R^ 
qS 


(10) 


Evaluation  of  the  coupling  function  C(wr,or)  requires  the  use  of  Gamma 
function  with  complex  argument.  A  series  expansion  [7]  given  by 


.  -z  z-1/2,  .l/2r,  1  1 

T(z)  =  e  z  (2ir)  [l  +  jj- +  - 


139 


571 


12z  288z2  51840z3  2488320z4 


163879 


5246819 


534703531 


209018880z5  75246796800z6  902961561600z 


7]  > 


(11) 


is  made  use  of  in  evaluating  the  hypergeometric  function  F. 
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Results  and  Discussion 

The  ratio  in  Eq.  (10)  which  does  not  contain  the  simplifications  of 
Ref.  [1]  is  evaluated  using  the  series  expansion  for  T(z).  The  result  is 
plotted  as  a  function  of  frequency  wr  in  Fig.  1.  The  amplitude  plot,  viz., 

20  log(h/hqg]  is  flat  up  to  the  break  frequency  and  then  falls  off  at 
approximately  3  dB  per  octave.  Up  to  about  the  break  frequency,  the  relation 
between  the  flow  rate  and  pressure  rise  can  be  expressed  by  the  usual  steady 
state  characteristic  curve.  On  the  other  hand,  at  frequencies  above  <dr  =  70, 
the  fluctuation  of  pressure  rise  is  negligible  even  though  the  flow  rate 
fluctuates  periodically.  An  observation  of  the  magnitude  plot  reveals  that  it 
is  similar  in  shape  to  that  of  a  first  order  control  lag  element  (of  the  form 
K/(xs  +  1))  with  a  break  frequency  at  oir  =  0.7.  However,  the  phase  plot  of 
the  dynamic  characteristics  (Fig.  1)  is  quite  different  from  that  of  a  first 
order  system.  The  phase  angle  (between  the  flow  fluctuation  and  pressure  rise 
fluctuation)  decreases  as  the  frequency  increases  until  the  break  frequency; 
then  it  increases  and  reaches  an  asymptotic  value  of  zero  for  large 
frequencies . 

In  the  proposed  transient  analysis  of  the  water  tunnel  circuit,  various 
components  of  the  circuit  that  influence  the  flow  dynamics  are  modeled  as 
point  elements.  As  noted  in  the  present  Introduction,  the  transfer  function 
for  each  element  is  then  utilized  in  writing  the  boundary  conditions  for  that 
element.  Since  the  dynamic  characteristics  of  the  tunnel  drive  pump  are 
available  from  the  present  analysis,  the  pump  can  be  fully  represented  in  such 


a  transient  analysis. 
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